Instability of higher dimensional charged black holes in the de Sitter world 
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We have shown that higher dimensional Reissner-Nordstrom-de Sitter black holes are gravitation- 
ally unstable for large values of the electric charge and cosmological constant in D > 7 space-time 
dimensions. We have found the shape of the slightly perturbed black hole at the threshold point 
of instability. Why only D = 4, 5 and 6 dimensional worlds are favorable as to the black stability 
remains unknown. 
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Introduction. The issue of stability of black holes 
was addressed for the first time yet in 1956, in a seminal 
paper of Regge and Wheeler [l| , who showed that four- 
dimensional Schwarzschild black holes, are stable against 
gravitational perturbations. This result confirmed that 
the Schwarzschild solution can indeed describe neutral, 
non-rotating black holes, because gravitationally unsta- 
ble systems simply could not exist. Later, the stability 
analysis was generalized for the Reissner-Nordstrom and 
Kerr solutions, which describes the electromagnetically 
charged non-rotating and neutral rotating black holes 
0- In 1992, the stability of asymptotically de Sitter 
black holes was proved by Mellor and Moss This 
meant that the general relativistic description of black 
holes is compatible with the idea of the expanding, de Sit- 
ter universe. As in four dimensional space-times, there is 
the uniqueness theorem for Kerr and Reissner-Nordstrom 
solutions, they were generally accepted as most physi- 
cally relevant. 

Last decade, the physical background has considerably 
changed with appearance of theories, implying existence 
of extra dimensions in nature, called brane- world theories 
@ , [H . These theories suggest a solution of the so-called 
hierarchy problem, that is the difference in scales of grav- 
itational and electro-weak interactions. In the scenario 
with Large Extra Dimensions 5], the (3-1- l)-dimensional 
brane, is embedded in a (4 -I- n)-dimensional space-time 
with n space-like compact dimensions. All matter is lo- 
calized on the (3 -I- l)-brane, while fields, which do not 
carry charge according to the Standard Model gauge 
group, can propagate in the bulk. An exciting oppor- 
tunity, that the brane-world theories give, is the possi- 
bility to observe the effects of strong, quantum gravity 
in a laboratory experiment at Tev energies. According 
to these theories, a miniature black holes might appear 
in the forthcoming experiments with particle collisions at 
the Large Hadron Collider or in the Cosmic Showers. 

When the black hole radius is much smaller than the 
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characteristic size of extra dimensions, one can describe 
the black hole by the Schwarzschild- Tangherlini metric 
When charged particles collide, a charged black holes 
must be formed. At the same time, recent observational 
data suggests the non-zero values of the cosmological con- 
stant in the Universe, so that the non-vanishing vacuum 
energy of the world must influence the formation of black 
holes. Thus, more general black hole background would 
be the Reissner-Nordstrom-de Sitter (RNdS) generaliza- 
tion of the Schwarzschild- Tangherlini metric. Yet, there 
is no traditional uniqueness theorem for D > 4-space- 
times, so that the important physical criteria that se- 
lects from all higher dimensional "black" objects (such 
as black holes, string, branes, rings, and saturns) is their 
stability: unstable objects cannot exist or need some 
mechanism of stabilization. 

Nevertheless, the stability analysis of D > 5 black 
holes became feasible relatively recently [8l|.~[l0l|. pdl ]. 
This reduction was performed for the D-dimensional 
Reissner-Nordstrom-de Sitter black holes in @] in the 
general form. Yet, the stability of the Reissner- 
Nordstrom black holes was proven analytically only for 
Z3 = 4, 5 space-time dimensions [Sj. The perturbation 
equations can be treated separately for all three types, 
called scalar, vector and tensor, according to the rota- 
tion group on the {D — 2)-sphere. When Z? — 4, we 
know the scalar type as polar and the vector type as 
axial, while the tensor type is usually a pure gauge. 
The higher dimensional cases were addressed in our ear- 
lier paper [l3|, where the stability of the D-dimensional 
Schwarzschild-de Sitter black holes was proved. Recently 
the stability of Reissner-Nordstrom-anti-de Sitter black 
holes (without dilaton) was shown in [l^. In addition, 
in 12] the numerical data for the quasinormal modes 
for vector and tensor types of gravitational perturba- 
tions of Reissner-Nordstrom-de Sitter (RNdS) black holes 
was given. Yet there it was claimed erroneously that 
Reissner-Nordstrom-de Sitter black holes are stable for 
all values of charge and A-term. In fact, in ^1^ for one 
particular, and most cumbersome, type of gravitational 
perturbations, the scalar type, one considered the effec- 
tive potential, which corresponds to the perturbations of 
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FIG. 1: The efltective potentials V- for p = 0.8, q = 0.9, 
£ = 2,3,4: (blue, green red respectively). As £ grows the peak 
becomes higher and the negative gap decreases. 



the Einstein equations with the frozen Maxwell field (see 
Eq. 8 in [13] )■ This approximation is valid when the 
charge of the black holes Q is considerably less than the 
black holes mass M, yet it is inappropriate for highly 
charged black holes. In the present paper we consider 
the dynamic behavior of the wave equation, which cor- 
responds to the complete perturbations of the Einstein- 
Maxwell equations, given by Eq. (5.61), (5.63 b) in 9]. 

Basic formulae. The metric of the D = d + 2- 
dimcnsional RNdS black holes is given by the line ele- 
ment 

ds^ ^ f{r)df ^ f-\r)dr^ ^r^d^d. (1) 

where dfl^ is the line element on a unit c?-sphere, /(r) — 
1~X + Z ~Y, 



X = 



2M 



Y = 



2Ar 



d{d+l)' 



Z = 



,,2d-2 ' 



(2) 



The equation of motion for gravitational perturbations 
of scalar type can be reduced to the wave-like equation 
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+ V± *(t,n) =0, 



(3) 



where the tortoise coordinate is defined as dr^ = 
dr/f{r) and the effective potential is a function of black 
hole parameters, r and of a multipole number £ that 
comes from separation of angular variables 



V±^V±{r,M,Q,K,D,£). 



(4) 



The explicit form of V± can be found in Q, formulas 
(5.61-5.63). The scalar type of gravitational perturba- 
tions, corresponding to the potential is the only type 
for which the stability cannot be proved analytically 
0. The potential V- reduces to the potential for pure 
gravitational perturbations, when Q — Q. On the con- 
trary, V+ reduces to pure electromagnetic perturbations 
propagating on the D-dimensional Schwarzschild back- 
ground in the limit of vanishing charge. 

We shall imply that 5* ~ e~*"*, lo = ujne — iujim, so 
that ujim > corresponds to a stable (decayed) mode, 



while U!im < corresponds to an unstable (growing) 
mode. If the effective potential V{r) is positive definite 
everywhere outside the black hole event horizon, the dif- 



ferential operator d /dri 



is a positive self-adjoint 



operator in the Hilbert space of the square integrable 
functions of r* , and, any solution of the wave equation 
with compact support is bounded, what implies stability. 
An important feature of the gravitational perturbations 
is that the effective potential F_ (see Fig. [1]), which 
governs the scalar type of the perturbations, has nega- 
tive gap for the lower values of the multi-pole numbers 
£. Higher i simply increase the top of the potential bar- 
rier, and are usually more stable lj|. Thus, we shall 



check here those values of £, for which the negative gap 
is present, and therefore the stability is not guaranteed. 

Numerical Method. We shall study the evolution of 
the black hole perturbations of scalar "-" type in time do- 
main using a numerical characteristic integration method 
that uses the light-cone variables m = t — and 
V — t ~\- r^. In the characteristic initial value problem, 
initial data are specified on the two null surfaces u = uq 
and V = vq. The discretization scheme we used, is 



^'(iV) = *(T4^) + *(£;)- ^'(S') - 

^, ViW)^{W) + ViE)^iE) 



(5) 



where we have used the following definitions for the 
points: N = {u+A, v+A), W = {u+A, v), E = {u, v+A) 
and 5* = {u,v). This method was very well tested 
for finding accurate values of the damped quasinormal 
modes (see for instance [l^ and references therein). Re- 
cently it has been also adopted for finding unstable, grow- 
ing, quasinormal modes in \yf\ for black strings, and in 
^14] for Gauss-Bonnet black holes. The agreement be- 
tween the time domain and accurate Frobenius methods 
is excellent. To test the reliability of the method, we in- 
creased the precision of the whole numerical procedure 
and decreased the grid of integration: unchanging of the 
obtained profiles of 4" signifies that we have reached suf- 
ficient accuracy of the computation. 

For convenience, we shall measure all quantities in 
units of the event horizon r+. Since the value of the 
event horizon is r+ = 1, the black hole mass is fixed as 
2M = 1 -I- Q2 _ 2A/d{d + 1). The cosmological constant 
written in terms of the cosmological horizon Tc allows to 
introduce the convenient variable p = r^jr^ = l/r-c < 1, 
so that 



2d(d+l) (1 + Q2)(/-1_1) 



„d+l 



1 



(6) 



We shall consider also the charge normalized by its ex- 
tremal quantity q — Q/Qext < 1- 

Discussion of the results. Let us start from the 
pure Reissner- Nordstrom (RN) black holes [p = 0). From 
the Table I one can see that quasinormal modes of non- 
extremal pure RN holes are damped for D = 6,7, ..11. 
For the near extremal values of charge Q, a power-law 
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TABLE I: I = 2 fundamental frequencies of gravitational perturbations of scalar type VL of D-dimensional RNdS black holes. 
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D = 


= 9 
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= 10 
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= 11 
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0.401i 
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1 
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0.863i 


2 AGO 


- 0.987i 


2.902 


- 1.087j 


0.1 
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254i 
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0.40H 
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608j 
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2.459 


- 0.983i 
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- 1.083i 
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0.400i 
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623j 
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0.822i 
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271i 
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0.412i 


1 


522 - 





596i 
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0.737i 


2.342 


-0.851i 
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0.8 
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260i 
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0.399i 
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490 - 





575j 
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252i 
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0.387j 
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459 - 
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855 - 


0.707i 


2.268 


- 0.823i 


2.689 


- 0.926j 


0.98 
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245i 
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0.380i 


1 
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552j 
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FIG. 2: Time-domain profile of near extremal q — 0.999 
Reissner-Nordstrom black hole perturbation {D = 11, p — 0). 
At the late time the power-law tail is observed (straight line 
in the logarithmic scale). The epoch of the quasinormal os- 
cillations becomes shorter for near extremal Q. 



FIG. 4: The parametric region of instability in the right upper 
corner of the square in the p — q "coordinates" for D — 7 (top, 
black), D = 8 (blue), D = 9 (green), D = 10 (red), D = 11 
(bottom, magenta). 
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FIG. 3; Time-domain profile of near extremal Reissner- 
Nordstrom black hole perturbation {D = 11, p — 0.8). q=0.4 
(brown) q=0.5 (blue) q=0.6 (green) q=0.7 (orange) q=0.8 
(red) q=0.9 (magenta). The smaller q, the slower growth of 
the profile is. 

damped tail dominates at asymptotically late times (Fig. 
[2]). When approaching near extremal Q, the epoch of 
quasinormal oscillations becomes much shorter (Fig. [2]), 
so that it is difficult to find out the accurate values of the 
QN frequency from the time domain profile, especially for 
higher D. Therefore some values in Table I are absent. 

The RNdS black holes are characterized by non-zero 
values of p and q. After careful testing of all the range 
of values of parameters q and p, we have found that for 



sufficiently large values of the both parameters, D > 7 
RNdS black holes are unstable. The typical picture of 
developing of instability can be found on Fig. [3l There 
one can see that for moderate values of q, the profile 
consists of damped quasinormal oscillations. Then, as 
the charge increases, the real oscillation frequency of the 
ringing decreases, approaching zero in the threshold point 
of instability. This is well understood, because unstable 
modes must be pure imaginary and the threshold point 
of instability corresponds to some static solution uj = 
of the wave equation [iTj . This is a natural picture for 
instability developed at lowest multi-poles. Instability in- 
duced by large i, on the contrary, appears as the growing 
of ^' after a long period of damped oscillations 'M\ . 

The parametric region of instability is shown on Fig. 
m For larger D the region of instability is bigger. An- 
other interesting question, which was beyond the scope of 
our paper, is if the extremal D = 6 Reissner-Nordstrom- 
de Sitter black holes are stable? Within the numerical 
method we can approach quite near the extremal values, 
but not the exact extremal limit. 

Black hole deformation at the edge of stability 
Let us study in detail situation when the charge is near 
the threshold of stability. We have found that there is a 
static mode (uj = 0) for £ = 2 "-" type of perturbations, 
while all other oscillations decay. We shall find here the 
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form of the perturbed metric at the threshold point of 
instabiUty in the hnear approximation as a background 
metric plus small perturbation, which is approximated 
by the found static mode. The correction to the black 
hole metric is given by Eq. (5.6a) of 9] 

Sgab = /ah(r-)S(z), Sgai = r/a(r)Si(z), 



Sg,, = 2r\HL{rH,§{z) + ffT(r)§.y (z)), (7) 

where S(z) are the scalar harmonics, which satisfy 

A§(z) + l)§(z). (8) 

A unit d-sphere is defined by the metric 

(if2^ = "fij{z)dz^dz\ 

with the angular variables {z}. 

From the Eq. (5.37) [| we see that = when uj — 
0. Consequently, in the Regge- Wilier gauge {fa{r) — 0, 
Hrir) = 0) we can see that the perturbed metric is static 
because /oi(r) = (see Eqs. (5.7b) and (5.8) in f^). 
Thus, the event horizon surface is given by the equation 



(R, If) plane 



5oo(0 + /oo(r-)§(2) =0. 



(9) 



Expanding this equation near the unperturbed black hole 
horizon r+ we find the expression for the black hole hori- 
zon f+ of the perturbed black hole 



5oo(^+) 



The radial coordinate R of the perturbed black hole is 
related to the coordinate r as 

i? = rVl + 2HL{r)Siz) w r(l + /7L(r)S(z)). (11) 

Substituting (fTO|l into (fTTjl we find for the perturbed hori- 
zon — + A3{z), where A is proportional to the 
amplitude of the linear perturbation 

A^foo{r+)/f'{r+)+r+HL{r+). 

We see that the null hypersurface is the sphere cor- 
rected by a small perturbation proportional to Si^2iz)- 
This surface can be easily imagined, because we al- 
ways can choose such an angular coordinate </> (associ- 
ated with the largest azimuthal number m = 2) that 
^e=2{z) = S^=2(0) and Eie=2{(t>) oc cos(2(/)) + const. Fig- 
ure [5] of the shape of the perturbed black hole at the 
instability suggests that the black hole could either split 
into two black holes or to transform into another solu- 
tion, such as, for instance, a black ring. The splitting 
into two black holes looks nevertheless improbable, be- 
cause such smaller holes would be themselves unstable 
and overcharged. 

One should note that since we consider a linear ap- 
proximation, we cannot describe the instability growth 




FIG. 5: The equatorial plane of the black hole horizon hyper- 
surfaces at the edge of stability. The dashed line corresponds 
to the unperturbed black hole of unit horizon. The blue, 
green and red lines correspond to the perturbed black holes 
after decay of all the dynamical modes. 



at later times because initially small perturbations will 
become comparable with the background at sufficiently 
late times. The final stage of the evolution of a unstable 
state must be described by the fully non-linear theory. 
The linear theory provides us only an initial guess for 
the developing of instability. 

Conclusions Let us enumerate the obtained results. 

1) The Reissner- Nordstrom black holes are stable for 
D = 5, 6, ..11. 

2) The D > 7 Reissner-Nordstrom-de Sitter black holes 
are unstable if values of the black hole charge and mass 
are large enough. 

3) The threshold values of parameters q and p, for 
which the instability appears, correspond to the domi- 
nance of some static solution of the wave equation. 

4) The shape of the slightly perturbed black hole at 
the instability point has been found. 

Strictly speaking the above instability does not mean 
that the corresponding black holes cannot exist, for in- 
stance a universe with considerable values of the cosmo- 
logical constant (presumably our universe in the early 
epochs), should have the chromo-dynamics instead of 
the considered here U(l) electrodynamics. If neverthe- 
less one wants to preserve U{1) electrodynamic, higher 
dimensional asymptotically de Sitter world, and large 
number of extra-dimensions D > 7 at the same time, 
he should take into account the above instability of black 
holes. This is important when investigating quasinormal 
modes or Hawking radiation of the Standard Model fields 
on the higher dimensional Reissner-Nordstrom-de Sitter 
background [l^l- The absence of coincidence of grav- 
itational and thermodynamic instabilities is also easily 
understood, if one notices that the thermodynamic insta- 
bility is induced by quantum effects, while gravitational 
instability is a classical effect. 
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